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STATISTICS  OF  ONE-DIMENSIONAL  CLUSTER  MOTION* 


John  D.  Wrigley,  David  A.  Reed,  and  Gert  Ehrlich 
Coordinated  Science  Laboratory,^  Department  of  Metallurgy,  and 
Department  of  Physics,  University  of  Illinois  at 
Urbana -Champaign,  Illinois  61801 


The  statistics  of  clusters,  made  up  of  metal  atoms  In 
adjacent  one -dimensional  diffusion  channels,  are  developed 
quantitatively.  Kolmogorov's  equation  Is  used  to  find  the 
mean  square  displacement  for  clusters  capable  of  existing  In 
energetically  different  configurations  at  the  same  displace- 
ment of  the  center  of  mass;  this  Is  done  under  steady  state 
conditions,  for  which  the  probability  of  finding  a specified 
configuration  does  not  vary  In  time.  Two  systems  are  examined; 

(1)  Dimers  capable  of  existing  In  an  Infinite  number  of  states, 
a situation  realized  it  dissociation  is  allowed,  and  (2)  Trlmers 
diffusing  on  planes,  such  as  W(211),  on  which  nine  distinct  jump 
processes  may  contribute  to  the  cluster  motion.  In  dimer  diffusion. 
It  is  demonstrated  that  dissociation  may  be  Important  even  If 
the  fraction  dissociated  Is  minor.  For  trlmers,  previous 
attempts  to  approximate  the  motion  through  the  use  of  average 
transition  rates  are  compared  with  the  exact  solutions  and 
found  wanting.  Important  statistical  quantities  beyond  the 


mean  square  displacement  are  presented  for  simple  dimers, 
capable  of  existing  In  only  two  states.  The  generating  func- 
tion Is  derived,  together  with  the  higher  moments  of  the 
displacement.  Probability  density  functions  for  the  number  of 
jumps  in  an  Interval  t,  for  the  waiting  time  up  to  a specified 
Jump,  as  well  as  for  the  displacements,  are  all  presented. 
These  differ  significantly  from  the  density  functions  for  an 
ordinary  random  walk.  However,  an  averaging  technique  allows 
simple  approximations  for  the  behavior  of  dimers. 
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Quantitative  observations  of  atomic  clusters  in  the  field  ion 

microscope  have  recently  made  it  possible  to  unravel  the  individual  jump 

processes  in  the  one-dimensional  surface  motion  of  dimers.^  This  analysis 

2 

rests  on  a stochastic  formalism  for  the  diffusion  process,  which  relates 

the  mean  square  displacement  of  the  center  of  mass  to  the  jump  rates  of 

3 

the  individual  atoms  in  a cluster.  Two  simple  statistical  problems 
have  so  far  been  considered  in  detail:  the  mean  square  displacement  of 

dimers,  made  up  of  atoms  in  adjacent  diffusion  channels  and  capable  of 
existing  in  the  two  states  labelled  0 and  1 in  Fig.  1,  and  the  diffusion 
of  model  trimers  for  which  three  jump  rates  suffice  to  describe  the 
motion.  Although  the  basis  for  an  atomic  analysis  of  cluster  motion  has 
been  established,  it  still  remains  to  elaborate  the  stochastic  formalism 
and  to  cover  the  behavior  of  more  complicated  real  systems.  That  is  the 
present  aim.  We  address  ourselves  to  three  related  areas: 

1.  In  the  diffusion  of  dimers,  dissociation  into  free  atoms  occasionally 

3 

occurs  at  high  temperatures.  To  assess  the  significance  of  these  events 
in  the  analysis  of  dimer  diffusion,  the  mean  square  displacement  is 
evaluated  for  dimers  that  can  exist  in  an  infinity  of  different  states. 

2.  The  actual  motion  of  trimers  may  depend  upon  nine  different 

4 

atomic  steps.  So  far,  this  complicated  motion  has  only  been  approximated. 
Here  the  diffusion  of  trimers  will  be  examined  in  its  generality,  as  a 
preliminary  to  the  analysis  of  experimental  observations,  allowing  all 
configurations  chat  can  be  achieved  by  moving  an  atom  one  lattice  spacing 
away  from  its  neighbor. 


2 


3.  Even  for  the  motion  of  the  simplest  cluster,  s non-dlssoclatlng 
dimer,  only  the  mean  square  displacement  has  been  worked  out.  For  a 
proper  understanding  of  cluster  diffusion,  as  well  as  for  assessing  edge 
effects  In  real  systems,^  the  distribution  governing  the  displacements  of 
dimers  must  also  be  known.  In  the  last  section  this  distribution, 
together  with  other  important  statistical  functions  defining  the  random 
motion  of  dimers,  is  developed. 


I.  FORMALISM  FOR  CLUSTER  DIFFUSION 

2 

In  the  previous  discussion  of  cluster  motion  in  one  dimension, 
only  systems  with  a unique  cluster  configuration  at  each  position 
of  the  center  of  mass  were  considered.  This  restriction  la  not  Inherent 
In  the  formalism,  however,  and  can  readily  be  removed  to  more  adequately 
describe  the  actual  surface  motion  of  clusters.  Throughout  this  presen- 
tation we  still  confine  ourselves  to  clusters  composed  of  atoms  In 
adjacent  diffusion  channels,  and  moving  In  one  dimension  only,  on  a surface 
with  spacing  1. 

We  proceed,  just  as  In  the  past,  to  follow  the  migration  of  an 
m-mer  by  mapping  the  location  of  Its  center  of  mass  on  an  Infinite  one- 
dlmenslonal  lattice,  with  spacing  f/m  and  m sites  per  unit  cell. 

The  random  walk  always  begins  In  the  unit  cell  based  on  x>0.  Diffusion 
of  this  cluster  Is  assumed  to  occur  one  atom  jump  at  a time;  regardless 
of  the  cluster  configuration,  each  jump  Is  presumed  to  span  only  a 
single  surface  spacing.  In  specifying  the  probability  of  finding  a 
point  X on  the  center  of  masa  lattice  occupied.  It  la  now  necessary  to 
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explicitly  recognize  the  existence  of  energetically  different  cluster 
configurations  jL,  The  probability  p that  point  x be  occupied  Is  now 


p ■ Z p . 
X ^ "^x.l 


1 - 0,1,2,. 


(1) 


where  the  summation  extends  over  all  configurations  ^ accessible  at  x. 

The  transition  rates  from  one  site  to  the  next  also  depend  upon  the 
configurations  Involved.  A general  formalism  to  handle  this  Is  outlined 
In  the  Appendix.  Here  It  suffices  to  specify  the  rates  starting  from  a given 
configuration  t,  with  X . denoting  the  rate  of  transitions  from  x to  x+1, 

X , 1 

and  M-  . the  rate  from  x to  x-1.  Kolmogorov's  equation^  for  the  probability 
p therefore  becomes 


dZ  p 


x,l 


dt 


“ ^ ^x-l,lPx-l,l  " ^ ^^x,l  ^x.l^^x.l  ^ ^x+l.l^x+l.l  ' 


with  the  additional  definitions 


X p = Z X .p  . 
x*^x  ^ x,l“^x,l 


^Px^^^,lPx,l  ' 


equation  (2)  reduces  to  the  usual  form 


“ip. 


~dT  “ ^-lPx-1  • ^^+%>Px  + '"x+lPx+l  “ 0«±l.±2.-.« 


(4) 
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The  differential  eqtiatlon  for  ((Ax)  ),  the  variance  of  x,  follows 
Just  as  previously 


- 2[<xX  > - <xu  > - <x)«X  ) - <u  ))]  + «X  ) + <u  » 

QV  XX  XX  XX 


(5) 


» 
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However,  the  possibility  of  having  different  cluster  configurations,  each 
with  its  own  characteristic  transition  probabilities,  at  any  point  x,  is 
now  explicitly  recognized  by  noting  that 


<X  > - E E X ,p  , 
X ^ ^ x,i'^x,i 


<M.  ) * E E M,  .p  . 

X . x,ix,i 


(6) 


X 1 


(xX)“ExEX  ,p  , 

X . x,i'^x,i 

X i ’ 


<X4^> 

X 1 


( x)  •»  E X E p 

X 1 


X,1 


Although  it  appears  that  examination  of  cluster  diffusion  in  terms 
of  individual  transition  probabilities  entails  greater  algebraic  com- 
plexity when  energetically  different  clusters  are  possible  at  a given 
site,  the  conceptual  framework  remains  essentially  the  same  as  that 
already  useful  in  earlier,  simpler  analyses.  As  specific  examples  of 
the  application  of  Eq.  (5)  to  cluster  motion,  we  will  analyze  the  linear 
diffusion  of  dimers  as  well  as  trimers.  In  experimental  studies,  the 
emphasis  has  been  upon  measurements  under  steady  state  conditions,  such 
that  the  probability  of  observing  different  cluster  configurations 
remains  unchanging  in  time.  In  view  of  the  complexity  of  the  systems 
of  interest,  we  shall  throughout  the  first  three  sections  limit 
ourselves  to  this  steady  state. 


II.  DIFFUSION  OF  DISSOCIATING  DIMERS 


A.  Mean  Square  Displacements 

In  considering  the  coupled  motion  of  two  atoms  in  adjacent  diffusion 
channels,  it  is  expedient  to  label  each  energetically  distinct  dimer 
configuration  by  the  separation  of  the  two  atoms,  expressed  in  units  of 
the  surface  spacing  projected  along  the  diffusion  direction,  as  in  Fig.  2, 
On  an  infinite  lattice,  an  infinity  of  configurations  is  therefore  possible 
on  each  site  x.  The  finite  range  of  interatomic  forces  insures  that  only 
a limited  set  of  these  configurations  corresponds  to  bound  pairs;  the 
rest  are  independently  migrating  adatoms.  This  distinction  will  subse- 
quently be  of  interest,  but  is  not  necessary  at  the  moment  - all 
configurations  can  be  treated  equivalently. 

The  rates  Important  for  dimer  motion  are  designated  according  to  the 
conventions  illustrated  in  Fig.  2.  Transitions  for  configuration  ^ to  i+1 
are  denoted  by  transitions  from  higher  to  lower  configurations,  that 
is,  from  jL  to  ^-1,  by  3^^.  The  scheme  of  rate  constants  connecting  various 
dimer  configurations  is  shown  in  Fig.  3.  Three  features  are  important; 

1)  Even  configurations  appear  on  even  sites  of  the  center  of  mass 
lattice,  odd  configurations  on  odd  sites.  2)  Inasmuch  as  only  single 
jumps  are  ever  allowed,  at  most  two  rate  constants  enter  into  the 
movement  of  a given  dimer  configuration  from  one  site  to  a neighbor. 

Thus 


x,i  '^x. 


where  6^^^  is  Kronecker's  symbol.  3)  The  rate  constants  depend  only  upon 
the  dimer  configuration,  and  not  explicitly  upon  the  position  x. 


It  is  now  an  easy  matter  to  evaluate  the  different  contributions 


I 

I 

I 

I 

I 

I 

f 


to  Eq.  (5).  The  rate  constants  averaged  over  position  become 


(X^“Z(X  p n+X  ,p  1+X  _p 

X x,0  x,0  x,l*^x,l  x,2  x,2  ’ 


Here  we  have  adopted  the  notation 


’^i  " ^ Px  i ^ “ 0.+1.+2,. . . 

X ’ 

to  denote  the  probability  of  finding  a specified  configuration  1^,  regard- 
less of  position  on  the  center  of  mass  lattice.  The  other  terms  appearing 
in  Eq.  (5)  now  follow  immediately. 


(xX^> 


<H^>  - S + 3^(l-6^o)]Pj  - <X^> 

^*^^x,0^x,0  ^ ^x,l^x,l  ^ ^x,2**x,2 

+ (Qf24e2)<x>2 

S [c^j-te^(l-6^Q)]<x)^ 


+ . 


+ 


^ "'’x.k 

X 


(8) 


(9) 


(10) 


(11) 


(12) 
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denotes  the  space  average  of  x for  a specified  configuration  k;  depending 
upon  the  parity  of  the  configuration,  this  average  will  include  either 
even  or  odd  sites  exclusively.  Finally, 

Only  two  terms  contribute  to  the  differential  equation  for  the 
variance  of  x,  which  now  appears  as 


d( (Ax)  ) , 


We  shall  be  concerned  with  the  behavior  of  dimers  over  time  intervals 

long  enough  for  a steady  state  to  be  established;  the  terms  P - E p 

i . "^x,i 

then  no  longer  vary  in  time.  Under  these  circumstances  the  requirements 
of  detailed  balance  are  satisfied  if 


!l  “O  '“2  “l  ‘■l+l  “l 


and 


’'O  »1  ’ >“1  ®2  ’ ■'l  » 


“l  “cPl 

" Q Pi  “ o D Pf>,  P, 


i+1 


i or : 


l -2  -1  V ■•t+1-’’o  " B- 


j»0  "j+1 


The  normalization  requirement 


( ! 


( 


ifi 


now  allows  us  to  write 


E P - 1 
i ^ 


p.  - (1  +E  n 


i J-0 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 
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In  Che  steady  state,  only  constants  appear  on  the  right  hand  In  Eq.  (14). 
Furthermore,  for  the  long  time  intervals  characteristic  of  the  steady 
state  the  precise  starting  point  becomes  unimportant.  Equation  (14) 
can  now  be  integrated  to  give  the  variance  of  x: 


<(Ax)  ) 


2t  E [a^-f8^(l-6.Q)]P^. 


Using  relations  (16)  and  (18),  this  becomes 


< (Ax)^) 


2t  p^q[1  + (a^-t83^)/3,^  + + . . .J 


4t  a. 


(1  + E n a,/e.) 

i j»l  ^ ^ 

1 

(1+E  na/g  ) 

i j»o  ^ J 


We  now  have  quite  a general  relation  for  the  variance  of  the 
displacements  executed  by  a dimer.  This  relation  is  capable  of  describing 
contributions  from  an  unlimited  number  of  different  configurations, 
provided  only  that  the  transition  rates  and  are  all  known.  For 
our  purposes,  however,  it  is  convenient  to  distinguish,  somewhat 
artificially,  between  two  types  of  configurations:  those  with  i small, 

for  which  Che  two  atoms  of  Che  dimer  are  bound  to  each  other,  and  those 
with  ^ large,  for  which  atomic  interactions  have  decayed  Co  such  an 
extent  that  the  atoms  move  independently.  Assume  that  the  latter  condition 
holds  for  1.  ^ f.;  then 


a 


1 


a 


(19) 


(20) 


1 a f. 


(21) 


where  ot  denotes  half  the  total  jump  rate  of  a single  atom.  The  factor 
of  1/2  accounts  for  the  fact  that  In  the  diffusion  of  single  atoms,  the 
total  jump  rate,  rather  than  the  rate  in  a specified  direction,  is 
traditionally  specified.  Equation  (20)  can  now  be  written  in  the 
alternative  form 


f-1 

2t{  E [^1  +e^(l-6l0)]P,  +2oP^}, 
i«0 


where 


P 

s 


E P 
1 s f 


i 


denotes  the  probability  of  finding  a free  dimer. 


(22) 


(23) 


B.  Specific  Examples 

In  the  limit  Pg*l,  that  is  as  dissociated  dimers  dominate,  the 
variance  of  the  displacements,  Eq.  (22),  becomes 

2 

( (Ax)  > ■ ^ort  . (24) 

This  describes  the  movement  of  two  free  atoms  on  the  center  of  mass 

lattice,  with  a grid  spacing  of  ^/2.  For  free  atoms,  an  expression  for 
2 

((Ax)  ) can  be  obtained  independently  of  the  preceding.  The  coordinate  x 
of  the  center  of  mass  is  related  to  the  coordinates  y^^  and  y^  of  the 
2 atoms  on  the  real  lattice  by 


X 


yi  + y-j 


(25) 


It  therefore  follows  that 


10 

If  the  movements  of  the  atoms  are  uncorrelated,  then 

<(Ax)^  - 2((Ay)^>  . 

For  the  random  walk  of  a single  atom,^  however,  we  have 

( (Ay)^>  ■ 2 a t , 

2 

yielding  for  ((Ax)  ),  the  variance  of  the  center  of  mass  coordinate, 
exactly  the  same  expression  as  derived  from  our  general  formalism  for 
pairs. 

If  dissociation  is  unimportant,  so  that  Pg“0  and  only  the  two  bound 
states  0 and  1 are  possible,  then  and  Eq.  (22)  transforms  to 

<(Ax)^>  - Aof^^t/(aQ  + 3j)  , 

the  relation  already  derived  in  I,  This  model  appears  to  describe  the 

available  data  for  dimers  on  the  (211)  plane  of  tungsten.  However,  in 

3 

quantitative  studies  by  Stolt  et  al.,  the  dimers  were  observed  to 
dissociate  occasionally.  In  order  to  estimate  the  effects  of  such  events, 
we  assume  that  dimers  in  configuration  1^3  are  entirely  free  and  that 
atomic  interactions  have  ceased  to  be  Important  once  the  interatomic 
separation  amounts  to  2Jt . This  situation  is  indicated  schematically  by 
the  potential  diagram  in  Fig.  4.  The  transition  rate  Is  then  just 
equal  to  the  rate  a in  the  free  dimer, 

or  - exp  - AE^/kT, 

where  AE^  is  the  activation  energy  and  the  prefactor  for  jumping. 


(27) 


(28) 


(29) 


(30) 


11 


The  variance  of  x now  becomes 

<(Ax)^>  - 4tPj^{Sj  +a^  +a[P2/(2p^)  + Pg/P^J]  . (31) 

The  probability  ratios  P9/P,  snd  P /P  depend  only  upon  the 
thermodynamics  of  the  different  dimer  configurations  and  in  principle 
are  accessible  to  direct  observation.  The  rate  a is  available  from  separate 
measurements  on  the  diffusion  of  single  atoms.  This  still  leaves  the  two 
transition  rates  3^^  and  with  only  Eq.  (31)  to  connect  them  to 
experiment.  However,  the  rate  at  which  dimers  change  from  configura- 
tion 1 to  2 is  fixed  within  limits  by  the  ratio  ^^e 

thermodynamics  of  configuration  2 relative  to  1.  In  the  worst  case,  we 
are  likely  to  have 

0^,  - V exp  - Ae"*"  /kT,  (32) 

“1  “1 

with  V ~ V , AE^  ■ AE^ ; a more  likely  event  is  that 

Of  O'  ^ 

AE*  ~ Ae"*"  + cE„  0 < c s 1,  (33) 

or  2 > V ' 

where  E2  stands  for  the  energy  of  configuration  2 relative  to  that  of 
configuration  1, 

Numerical  estimates  for  the  effects  of  dissociation  upon  the  mean 

square  displacement  are  shown  in  Fig.  5 for  a system  modeled  on  the 

3 

behavior  of  rhenium  dimers  on  W(211).  It  appears  that  dissociation  can 
significantly  affect  the  mean  square  displacements.  If  a sizeable 
fraction  of  the  dimers  can  transgress  beyond  the  0 and  1 configurations. 
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2 * 

then  ((Ax)  ) diminishes,  as  the  mobility  of  free  atoms  in  this  model 
system  is  less  than  that  of  bound  dimers.  It  is  obvious  that  these 
effects  must  be  allowed  for  in  the  analysis  of  experimental  observations. 

The  corrections  for  dissociation  do  not  sensitively  depend  upon  the 

particular  assumptions  concerning  the  unknown  rate  of  transition  from 

configuration  1 to  2.  Even  in  the  unlikely  worst  case,  the  variance  of  x 

differs  only  by  ss  10%  from  that  found  with  the  more  likely  smaller  values 

of  for  which  the  actual  magnitude  of  the  rate  has  hardly  any  effect 
2 

on  ((Ax)  ).  It  is  of  interest  that  the  dependence  of  the  variance  upon 
the  particular  assumptions  concerning  is  greatest  when  the  fraction  of 
dimers  in  state  2 is  least.  In  principle,  it  is  possible  to  have  an 
entirely  negligible  amount  of  dissociation  and  yet  have  important  contri- 
butions to  the  diffusion  from  configurations  such  as  2.  Physically  this 
happens  if  transitions  from  configuration  1 to  2 are  possible,  but  the 
return  rate  is  much  faster.  No  appreciable  concentration  can  build  up 
in  state  2;  every  time  a transition  into  this  configuration  occurs,  it  is 
immediately  followed  by  the  reverse  process.  However,  each  of  these 
events  means  a displacement  of  the  dimer  above  and  bevond  that  of  a 
dimer  with  only  two  allowed  configurations. 

It  is  important  to  consider  this  possibility  in  interpreting  the 
diffusion  of  dimers.  Even  if  no  dissociation  is  actually  observed  in 
an  experiment,  the  analysis  should  properly  be  carried  out  with  a model 
incorporating  the  effects  of  dissociating  dimers.  This  can  be  accomplished 
using  the  thermodynamic  information  to  guess  at  the  value  of  the  transi- 
tion rate  from  configuration  1 to  2,  and  estimating  the  Jump  rates  with 
and  without  allowing  dissociation. 
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III.  TRIMER  DIFFUSION 


The  linear  migration  of  trimers  can  be  handled  in  much  the  same 
way  as  that  of  dimers,  although  the  nomenclature  becomes  more  complicated. 
The  unit  cell  of  the  lattice  on  which  motion  of  the  center  of  mass  is 
mapped  now  contains  three  points;  as  shown  in  Fig.  6,  two  energetically 
different  trimer  configurations  are  possible  at  each  such  point.  The 
convention  found  useful  in  designating  dimer  configurations  in  Sec.  II 
no  longer  suffices  for  trimers.  We  therefore  retain  the  labelling  devised 
in  I.  In  this  scheme,  slanted  or  V shaped  configurations  are  distinguished 
from  the  others  by  adding  B or  A respectively  to  the  usual  configuration 

9 

index.  The  different  terms  entering  Eq.  (5)  for  the  variance  of  x can  now 
be  readily  evaluated.  Keeping  in  mind  that  trimer  configuration  0 occurs 
only  on  sites  of  type  0,  configuration  1 on  sites  of  type  1,  etc.,  it 
iirir.jd lately  follows  that 


OB 


+ (bj+bjj)  E^p^^^^  + bjj^  E^p^^^g  + (c^4Cjj)  E^p^^^A 


(34) 


‘^III  ^2Px,2B  ; 


symbols  such  as  E^,  for  example,  indicate  summation  over  x for  all  sites 
of  type  1.  The  averaged  transition  rate  (X^)  can  be  written  more  compactly 

as 


^^x^  ■ ^*l‘^II^  ^QA  ^‘^l'^II^^2A 


“lII^OB  ‘’lII^lB  ®III**2B 


(35) 


where,  in  conformity  with  the  conventions  for  dimers,  we  use  the 
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abbreviation 

^QA  ^ ^QPx.OA 

to  describe  the  probability  of  finding  a configuration,  such  as  OA  in  this 
Instance,  regardless  of  position  x.  We  recall  that  under  steady  state 
conditions  these  probabilities  are  connected  by  a number  of  simplifying 
relations : 


"l 

P, 


2A 

P ■*  P 

IB  *^2B 

’oA 

P * P 

IB  '^OA 

PqA  - [1  + 2aj/cj  + 2ajj/Cjjj  + 

Equation  (35)  therefore  simplifies  to 


-1 


<X^>  . (aj  + + (bj  + bjj  + cj  + Cjj)P 


lA 


Similarly, 


■ <*I  <■>1  >'II  + 'l  + 'll>fl* 


<xX^>  + + (b,*j,)p^_ 


lA 


* ’’lIlPx,lB  ^‘"l'*^II^Px.2A  ‘"lIlPx,2B" 
+ •iii(x>0B  + **111^’'^  IB  ‘^111^*^28  • 


(36) 


(37) 


(38) 


(39) 


(40) 
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Here  we  have  resorted  to  the  convention 


(x)^  =2iXP^,1A 


(41) 


to  indicate  the  average  of  x for  a specific  trimer  configuration,  in  this 
instance  lA.  The  last  remaining  term  in  Eq.  (5)  is 


(42) 


The  differential  equation  for  the  variance  of  x can  now  be  written  as 


(Ax)^> 

dt 


- 2{(bj-Cj+bjj-Cjj)(<x)^^-<x>2^)  2b^ 

+ Po^[2(aj+aj^)  + (b^+2b^^+2Cjj)aj/Cj]}  . 


(43) 


We  still  require  the  two  terms 


^*\a  " ^*^2A’  h ^*\b  " ^’‘^2B  • 

To  derive  these  we  need  the  probability  p of  having  a specified 
configuration  at  the  point  x.  This  can  be  obtained  from  Kolmogorov's 
equation,  as  shown  in  the  Appendix.  For  p^  as  an  example,  the  appropriate 
differential  expression  is 

*^^x.lA 

dt  ' *l’’x-l,0A  *IIlPx-l,0B  ■ ^‘’I'^II^l'^II^Px.lA 

**lPx+l,2A  **IIlPx+l,2B  * 

It  follows  in  the  usual  way  that 


ai«x>o^+Po^)  - 


In  much  the  sami.  fashion,  It  can  be  shown  that 


+ «l(<’‘)aA"^QA^  ■*■  ®III^^’‘^0b"^0B^  ■*■  ^III^^*Hb‘‘’^1B^  • 
Using  the  steady  state  relations  (37)  and  the  definitions  of  Eq.  (44), 


we  f ind 


_ . . (2bj4Cj4bjj4Cjj)X^  - bjjjXg  - 2(bj-Cj4bjj-Cjj)Po^a^/Cj  . 

The  differential  equation  for  follows  from  an  analogous  sequence  of 


steps  as 


-71  ■ - • '’lA  • - ■ll'foA- 

Only  the  steady  state,  in  which  the  time  dependence  becomes  unimportant, 
is  of  interest.  Equations  (48)  and  (49)  therefore  reduce  to 

(2bi  + Cj  + bjj+Cjj)X^  + bjjjXg  - - 2(bj-Cj.+bj.j.-Cjj.)PQ^8j/Cj 


with  the  solution 


*l\ 

V • 2P  — ^ 

^A  OA  CjA 


^ 


OA  CjA 


(51) 
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4 • 


III 


Substituting  these  results  in  Eq.  (43)  and  integrating,  the  variance  of 
the  center  of  mass  displacements  for  trimers  appears  as 


+ Cjd+ajj/a^)  + 0>j+2bjj+2Cjj)/2]aj/Cj  . 


(52) 


In  the  limit  as  a^^  ■»  = c^.^  •»  0,  trimer  migration  can  only  proceed 

through  A configurations.  This  nossibility  has  been  treated  exactly  in  I. 
For  such  a diffusion  path,  Eq.  (52)  simplifies  to 


((Ax)^> 


(2aj.-K;j)(2bj4Cj) 


(53) 


vidiich  is  precisely  the  relation  previously  derived.  If  the  transition  rates 
meet  the  condition  then  diffusion  must  occur  by  movement  of  the 

trimer  through  the  configurations  OB,  lA,  and  2A  only.  This  constitutes 
another  symmetrical  path,  for  which  Eq.  (52)  yields 


<(Ax)^> 


(2aiii+Cii)(2bi+Cii) 


(54) 


again  in  agreement  with  previous  work. 

The  complete  equation  for  trimer  diffusion  involves  nine  transition 
4 

rates.  These  are  not  all  accessible  to  direct  observation  in  any 

straightforward  way,  but  Eq.  (52)  makes  it  feasible  to  model  the  behavior 

of  trimers  for  a known  interatomic  force  law.  An  approximate  procedure  has 

2 

been  suggested  in  the  past.  In  this  the  motion  is  Idealized  as  occurring 
via  a symmetric  path  through  the  unit  cell;  average  rate  constants  are 
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then  used  to  account  for  the  presence  of  energetically  different  con- 
figurations at  a lattice  point.  Tests  with  model  potentials  reveal  that 
this  approximation  works  surprisingly  well  in  many  instances,  yielding  mean 
square  displacements  within  25%  of  the  true  values.  Under  these  circum- 
stances the  diffusion  parameters  deduced  from  the  usual  Arrhenius  plots 
are  essentially  without  error.  However,  for  some  combinations  of  physically 
reasonable  rate  constants  the  average  rate  approximation  fails  by  more 
than  an  order  of  magnitude.  Despite  its  complexity,  the  exact 
relation  (52)  for  the  variance  of  the  center  of  mass  positions  must 
therefore  be  used  in  modeling  the  diffusion  of  trimers, 

IV,  STATISTICS  OF  TWO -STATE  DIMERS 

The  emphasis  up  to  this  point  has  been  upon  the  mean  square  displacements 
executed  by  the  center  of  mass  of  a cluster  during  surface  diffusion. 

The  distribution  of  displacements,  as  well  as  higher  moments  of  the  distribu- 
tion, are  also  of  concern.  From  a practical  point  of  view,  we  note  that  one- 
dimensional  diffusion  is  observed  on  planes  of  finite  size,  so  that  corrections 
have  to  be  made  for  the  possible  effects  of  the  edges.  The  simple  relations 
available  for  this,  however,  are  based  on  the  assumption,  as  yet  not  tested, 
that  the  displacements  can  be  approximated  by  a Gaussian  distribution. 

More  than  that,  the  one-dimensional  motion  of  clusters  also  constitutes 
a novel  type  of  random  walk,  of  interest  in  its  own  right.  Here  we  there- 
fore derive  several  statistical  quantities  Important  in  describing  cluster 


if 


motion. 
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We  confine  ourselves  to  the  simplest  clusters,  to  dimers  that  can  only 
exist  in  the  two  configurations  0 and  1 indicated  schematically  in  Fig.  1, 
Transitions  from  x>»0  to  x“l  are  assumed  to  occur  at  the  rate  a,  from  x=l  to 
x-0  at  the  rate  b.  That  is,  in  the  general  scheme  indicated  in  Fig.  2, 

Q'g“£  other  rates  vanish.  To  ease  the  derivation  of  higher 

moments,  which  in  principle  are  directly  accessible  through  Kolmogorov's 
equation,  we  first  develop  the  probability  generating  function  G(z),  defined  by 

G(z)  £ E z^p^  X = 0,  + 1,  + 2,...  . (55) 

X 

The  distribution  governing  dimer  displacements  is  presented  separately. 

Also  of  interest  for  characterizing  the  random  walk  of  dimers  is  the 
probability  density  function  for  the  total  number  of  jumps,  and  for 
the  time  required  to  make  a fixed  number  of  jumps,  the  flrrr  of  which  plays 
a key  role  in  the  distribution  of  displacements.  In  the  last  section,  both 
are  derived  and  compared  with  the  behavior  of  simple  random  walks. 

A.  Generating  Function  and  Moments 
Starting  from  Kolmogorov's  relation,  Eq.  (4),  it  follows  quite 
generally  that 


- E z*>  ,p  , - E z*(X  + u )p  +Ez’‘x  ,p 
dt  x-r  x-1  X X ^x  x+1 


x+1  x+1 


For  dimers,  this  differential  equation  reduces  to 


G(  ) - (u  - 2)(aGQ  + bG^), 


where  G(z)  denotes  the  time  derivative  of  the  penerating  function  and  we 
have  introduced  the  definitions 


(56) 


(57) 


u 


G 

G 


0 

1 


z 

E 

E 


+ l/z 


P 

0 *^x 


1^  P 
1 X 


(58) 
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To  proceed  further,  the  sums  Gq  and  Gj^  over  even  and  odd  sites^  respectively, 
must  be  expressed  In  terms  of  the  dimer  transition  rates.  Again  from 
Kolomogorov' s relation,  we  find  the  coupled  differential  equations 


S ■ '’“S  - 

— siiCq  “ 2bGj^ 

Laplace  transformation  converts  these  into  the  form 

(s  + 2a)X{GQ]  - baC[G^}  = Gq[  t-O] 

- auXlGg}  + (s  + 2b)X[G^}  = G^[t-0]  ; 

s is  as  usual  the  argument  of  the  Laplace  transform,  and  the  values  of 
Gq  and  G^  at  time  t=0  are  denoted  by  Gq[ t=0]  and  G^[t=0]  respectively. 
For  excursions  starting  from  even  sites, 

G^°^[t-0]  = 1.  Gj°>[t=0]=  0, 


for  excursions  from  odd  sites, 

,(1) 


G^*'^[t«=0]  = 0,  G^^Lt^O]  - z; 


(59) 


(60) 


(61) 


(62) 


the  origin  of  the  dimers  is  indicated  explicitly  by  the  appropriate  super- 
script. For  the  moment  confining  ourselves  to  excursions  starting  at  an 
even  site,  we  find 


(0)i 
0 ^ 

(0): 


s + 2b 


(s  + q)(s  + r) 


au 

(8  + q)(8  + r) 


(63) 


where 
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q - (a  + b){  1 + [l  + ab(u^  - 4)/(a  + b)^]^] 
r s (a  b){l  - [l  + ab(u^  - 4)/(a  + b)^]^}  . 

The  inverse  transformation^ of  Eq.  (63)  yields  the  desired  quantities 


(64) 


,(0)  ^ (q  - 2b) 
'0 

,(0) 


exp(.  qt)  - exp(-  rt) 


(65) 


[exp(-  qt)  - exp(-  rt)] 


(q 


Substituting  these  in  Eq.  (57)  and  integrating  gives  the  generating  function 

.(0) 


G (z)  for  starts  at  even  sites  as 


G^°\z)  = 1 + [[q-b(u+2)][l-exp(-qt)]/q 

- [r-b(u+2)][l-exp(-rt)] /r}  . 


(66) 


A similar  sequence  of  steps  yields  the  generating  function  G^^\z)  for 
excursions  from  odd  sites: 

G^^\z)  = z + [[q-a(u+2)][l-exp(-qt)]/q 

- [r-a(u+2)] [l-exp(-rt)] /r]  . (67) 

In  the  limit  as  the  transition  rates  ^ and  b become  equal,  the  generating 
function  for  starts  at  x^  reduces  to 

G (z)  « z expa(u-2)t  a-b,  (68) 

appropriate  for  a simple  random  walk  over  discrete  sites  but  continuous 


in  time. 


r 


I 

I 

I 

i 

I 

I 

I 

I 

I 

I 
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th  n 

The  n moment  (x  > of  the  coordinate  of  the  center  of  mass  is  now 
easily  accessible,^  through  repeated  differentiation  of  the  generating 
function  according  to 

(z  ^ . (69) 

Z"1  X 

We  display  only  the  2nd  and  4th  moments;  the  odd  moments  vanish  for 
excursions  starting  from  even  sites. 


2v(0) 


(x2)<l> 


^ 2(a^)  tl-exp-2(a+b)t]} 


a+b 


^ • 2(fiS  [ 1 -exp-2 (a-fb)t]} 


(70) 


/ 4s (0)  Sab  r , . 2 _ (a-bl 

(x  )'  » [6abt  +3at  ' exp-2(a+b)t 

(a-+b)'^ 

. (5a^-5ab+2b^)t  (a-b)fa^-28ab+7b^>  ..  . . 


1 + — {6abt^-3bt  ^ exp-2(a+b)t 


(71) 


(a4b) 

V“  ■ “ / 

(5a^4ab-f«b^)t  (a-b)(13a^-16ab+7b^l  n. 

+ ^ 2 [l-exp-2(a4b)t]} 

^ ’ Sa(a-fb)^ 


In  the  limit  £“b,  the  moments  simplify  to 

<x2)(0>  - 2at 

- 2at  + 1 (72) 

<x^<0>  - 12aV  + 2at 

• 1 + 12a^t^  + 14at  . (73) 

These  are  precisely  the  values  for  an  ordinary  random  walk  In  continuous 
time . 


For  a Gaussian  walk  characterized  by  a total  jump  rate  2X  and 
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starting  at  the  origin. 


/ ^\  •»  / 2v2 


<x  > - 2).t 


(74) 


We  note  that  In  the  limit  of  long  times,  only  the  first  term  In  braces 
occurring  In  Eqs.  (70)  and  (71)  contributes.  In  this  limit,  therefore, 
the  ratio  of  the  4^^  to  the  2*^^  moment  for  dimer  displacements  Is  exactly 
that  for  a Gaussian  walk.  The  2*^*^  moment  for  dimers  Is,  of  course,  equal 
to  that  of  a Gaussian,  but  with  an  effective  jump  rate  X ■ 2ab/(a+b). 

In  short,  for  the  lower  moments  at  least,  a Gaussian  approximation  Is 
appropriate. 

B.  Probability  Densities 

1.  General  Approach 

In  principle  the  distribution  governing  cluster  displacements  can  be 
obtained  from  the  generating  function  G(z)  by  expansion  In  a ! aurent 
series.  We  adopt  a simpler,  physically  more  transparent  approach. 

For  the  center  of  mass  to  end  at  point  x after  N jumps,  having 
started  at  x^.  It  must  take  [N  + (x  - Xq)]/2  steps  to  the  right  and 
[n  - (x-Xq)]/2  to  the  left;  Jumps  to  the  right  and  to  the  left  are  assumed 
equally  probable.  The  conditional  probability  f{x|N}  of  ending  at  x,  given  a 
total  of  N jumps,  therefore  obeys  the  binomial  law^ 


f{x|N} 


N-Hc-x 


(75) 


note  that  f{xlN]  ■ 0 If  N and  (x-x^)  are  not  of  the  same  parity.  The 
joint  probability  f(x,N}  of  the  center  of  mass  being  at  x after  N jumps  Is 


f{x,N}  « p^(t)f[x|N} 
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(76) 


where  Pj^(t)  denotes  the  probability  that  during  an  interval  t a dimer  will 
make  a total  of  N jumps. 

To  obtain  the  desired  probability  p of  a dimer  ending  at  x,  we  just 

12 

Sum  over  all  jumps  N 

00 

P„  “ 2 p (t)f{x|N]  . (77) 

N-0  ^ 


It  now  remains  to  derive  the  probability  density  P^(t)  describing  the 

total  number  of  jumps  executed  by  the  center  of  mass  during  time  t. 

Closely  connected  to  the  number  of  jumps  during  such  an  interval  is  the 

time  required  to  make  a fixed  number,  k,  of  jumps.  Although  not 

directly  involved  in  arriving  at  p^,  the  time  distribution  is  of  interest 

in  characterizing  those  aspects  of  the  motion  of  dimers  distinct  from 

ordinary  random  walks.  The  probability  densities  for  N and  for  Tj^  are 

therefore  derived  separately  in  the  next  two  sections. 

2.  Probability  Density  for  the  Number  of  Jumps  N 

We  seek  the  probability  Pjj(t)  that,  as  the  dimer  moves  along  the  lattice, 

it  will  make  a total  of  N jumps  in  the  time  interval  t.  Let  X denote 

N 

th 

the  rate  of  making  the  (N+1)  jump.  Inasmuch  as  the  number  of  jumps  N 

13 

can  never  decrease,  Kolmogorov's  equation  for  Pj^(t)  becomes 


dt 


(78) 


dPj,(t) 


dt 


N > 0.  (79) 
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The  Laplace  transforms  of  Pj^(t)  are  therefore 


I r 


-7^ 


N 


Equations  (80)  and  (81)  Immediately  lead  to  the  general  relation 


N > 0.  (81) 


•C[PN(t)} 


-L_  V . 

" + ^0  j-0  ® S+1 


N > 0.  (82) 


If  excursions  start  at  an  even  site,  then  the  rate  of  jumping  for  dimers  is 


given  by 


X^.2a 

X^.2b 


N even 


N odd 


We  shall  only  present  solutions  for  this  situation,  but  note  that  for 
excursions  starting  at  an  odd  site,  a and  b are  Interchanged  throughout. 
With  the  rates  defined  by  (83), 


j.r_(0),j.)]  - (2a)^^^(2b)^^^ 

( s+2a ) 


N even 


N odd 


For  an  odd  number  of  jumps  the  inverse  transform  is  immediately  available, 


4 » 


and  gives 


Pn 


. (2a)^^'*'\^^^(2b)^^'^^^^  nV^^expL-  (a-^)t]  1 . C(a-b)t] 


[2(a-b)]‘ 


N odd  . 


The  solution  for  an  even  number  of  jumps  Is  obtained,  by  substituting 
(86)  in  Eq.  (79)  and  integrating,  as 


pr<« 


„4^(N+1)/2 


2(N/2).' 


giSPi-:  ja-H))tjrj.  rra-bu1 

[2(a-b)](N“l>/2  ^ \n-1)/2L(® 


* T 
; i 


N even  , (87) 

is  the  modified  spherical  Vessel  function  of  the  first  kind. 

As  the  rates  ^ and  b approach  each  other,  dimer  behavior  becomes 
independent  of  site,  and  the  total  number  of  jumps  should  therefore 
conform^  to  a Poisson  distribution.  The  limiting  form  of  Pjj(t)  when  a »b 
can  be  obtained  making  use  of  the  series 


PL  T z“ 

'2z  n+1^  “ l*3*5...(2n+l) 


and  the  identity 


1 ,oos 

l.'(2n+3)  2,'(2n+3)(2n+5)  ^ 


l*3*5...(2k-l)*(2*^k:)-  (2k)! 


Both  Eq.  (86)  and  (87)  reduce  to  the  proper  form 


o - (2at)  exp_-_2,at 
V ^ N! 


The  probability  density  for  the  number  of  Jumps  by  a dimer  is 
compared  in  Fig.  7 with  the  behavior  typical  of  an  ordinary  random  walk, 
for  which  the  jump  rates  are  independent  of  position,  and  a ■ b . As  expected, 
when  Che  Jump  rates  differ  significantly,  there  are  pronounced  variations  in  the 
probability  density  on  going  from  an  odd  to  an  even  number  of  Jumps. 


L 


In  Fig.  8,  the  average  of  the  probabilities  for  two  neighboring  jumps, 
trhen  the  rates  a and  b are  significantly  different,  is  compared  with  the 
same  average  for  a Poisson  distribution,  characteristic  of  an  ordinary 
random  walk.  The  agreement  is  quite  reasonable.  Averaging  over  adjacent 
values  of  N eliminates  the  large  excursions  in  the  probability,  revealing 
overall  trends  for  dimers  which  are  quite  similar  to  those  in  a simple 
random  walk. 

3,  Probability  Density  for  the  Waiting  Time  T^ 

Basic  to  the  development  in  this  section  has  been  the  assumption  that 
cluster  diffusion  can  be  described  as  occurring  via  atomic  jumps  random  in 
time,  with  transition  rates  depending  only  upon  the  state  of  the  system 
from  which  the  jump  takes  place.  That  is,  the  elapsed  time  t between 
jumps  ji  and  i,  + 1 has  an  exponential  probability  density 

- Xj^exp  - X^T  , (91) 

where  X^  denotes  the  rate  of  the  (i+1)^*^  jump.  The  time  T^^  at  which  the 
k^^  jump  takes  place  is  given  by  the  sum 

k-1 

• Z J,  . (92) 

^ 1-0  ^ 

The  waiting  times  for  individual  jumps  are  Independent  quantities,  each 
with  a characteristic  function 

® X 

(j£(9)  -J  exp(jOr)f^(T)dT  - Tq  , (93) 

-CO  1 

2 

where  J - -1.  Because  of  this  Independence,  the  characteristic  function 

0 (0)  for  the  probability  density  of  T.  is^ 

T, 
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k-1  k-1  X 

(9)  - n (j  (0)  - n . (94) 

k 1-0  ^ 1-0  ^1  ■ 


Substituting  the  jump  rates  (83)  appropriate  to  dimers,  the  characteristic 
function  of  for  excursions  starting  at  an  even  site  Is  found  as 


4°^(9)  - [2a/(2a-j0)]  ^'‘■'■^^^^[2b/(2b-j9)] 
k 

k odd 

(95) 


(9) 


[2a/(2a-j8)]‘^^^[2b/(2b-je)]'^^^ 


k even. 


For  excursions  from  odd  sites,  a and  b must  be  Interchanged.  The  desired 

probability  density  function  ^ (t)  Is  now  obtained  from  the  exponential 
14  ^ 

Fourier  transform  of  0 (9), 

^k 


00 

k -»  k 


(96) 


as 

f^°\t)  - (2a)^*‘‘‘'^^^^(2b)  ^^/^exp(-2at) 

X [t^‘V(k-l)']^F^[(k-l)/2;k;(2a-2b)t] 


k odd 


(97) 

f (t)  - (2a)‘‘^^(2b)*'/2exp-(2at)[t‘‘‘V(k-l)']-F,[k/2;k;(2a-2b)t] 

k even  . 

Here  jFj^(a;  c;  z)  denotes  Kummer's  confluent  hypergeometrlc  series. 

In  the  limit  as  a-b,  jumping  becomes  a Poisson  process,  in  which  the 
waiting  time  for  the  k^^  jump  Is  known  to  conform  to  a gamma  distribution.^ 
Also  In  this  limit,  the  Kummer  functions  In  Eq.  (97) 
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approach  unity.  The  probability  denalty  for  the  time  at  which  the 
jump  occurs  then  Is 

(t)  - (2a)‘"[t'''V(k-l).']exp  - 2at  . (98) 

k 

As  expected,  this  is  just  the  gamma  distribution  of  order  k. 

Probability  densities  according  to  Eq.  (97)  are  plotted  in  Fig.  9. 

In  these  graphs,  the  probabilities  calculated  for  different  values  of  the 
rates  a and  b are  compared  with  each  other,  always  keeping  X » 2ab/(a+b) 
constant.  When  the  number  of  jumps  k is  small,  different  choices  of  £ 
and  ^ yield  significantly  different  curves.  Of  course  for  k even,  the 
plots  are  invariant  to  an  interchange  of  a and  Despite  that,  the 
probabilities  for  £ ^ are  not  similar  to  those  for  a * b.  As  k becomes 

large,  however,  these  differences  diminish;  the  waiting  times  can  then  be 
reasonably  approximated  by  a gamma  distribution,  calculated  for  the 
effective  Jump  rate  X. 

4.  Distribution  of  Displacements 

With  the  probability  density  Pj^(t)  for  the  number  of  jumps  during  an 
interval  t in  hand,  we  can  now  write  out  the  probability  density  for  x. 

Using  Eqs.  (75)-(77)  as  well  as  (86)  and  (87),  we  find  that  for  starts  at 
an  even  site. 


-~^^^|^^exp[  - (a+b)t]L 


N/2 

(4ab)  n; 


j.N/2 

^ 2^[2(a-b)  ‘ ^(N+l)/2'- ^ 


x,N  even 


(99) 


- /!k 


exp[-  (a-rt))t] 


N/2 

/ (4ab)  n; 


N/2 


"(¥K¥)<¥)' 


I.,„[(a-b)l;] 


x,N  odd 

If  the  initial  position  is  an  odd  site,  x on  the  rhs  of  Eq.  (99)  is  replaced 
by  x-1,  and  the  rate  constants  a and  b are  interchanged. 

In  the  limit  a»b,  it  follows  immediately  from  Eqs,  (88)  and  (99) 

that 

-1 


p^  = exp(-  2at 


) I (at)”  - 2at)I^(2at).  (100) 

N 


That  is,  the  probability  density  of  x reduces  to  the  form  known^  to  hold 
for  an  ordinary  random  walk  continuous  in  time.  It  is  of  interest  to 
compare  this  limiting  density  with  that  describing  the  actual  displacement 
of  dimers,  according  to  Eq.  (99).  This  is  done  in  Fig.  10.  The  distribu- 
tions clearly  show  two  effects:  1)  A gradual  decrease  in  the  probability  of 
being  at  a point  x as  x increases.  2)  Strong  oscillations,  caused  by 
differences  in  jump  rates  from  even  and  odd  sites.  It  is  clear  that  a 
Gaussian  can  be  seriously  in  error  as  an  approximation  to  any  individual 
value  of  p^.  Averaging  the  probability  over  adjacent  sites,  however, 
reduces  this  problem.  In  Fig.  11,  these  averaged  values  are  compared  with 
a Gaussian,  all  at  the  same  effective  rate  X 2ab/(a-ri5).  Even  when  the 
rates  a and  b are  quite  different,  a Gaussian  serves  as  a reasonable 
approximation  to  the  averaged  probability  of  finding  a dimer  in  a 


specified  unit  cell. 
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V.  SUMMARY 

Elementary  statistical  techniques  have  been  used  to  evaluate  exactly 
the  mean  square  displacement  of  small  clusters,  capable  of  existing  in 
energetically  different  configurations  at  each  site  of  the  center  of  mass 
lattice.  Two  specific  problems  have  been  examined  - the  diffusion  of 
dissociating  dimers,  and  of  bound  trimers.  The  techniques  used  for  the 
latter  should  also  be  adequate  to  describe  the  motion  of  larger  clusters. 
However,  the  present  development  already  suffices  for  the  analysis  of 
available  experimental  information,  which  will  be  presented  separately. 

The  evaluation  of  the  probability  densities  described  here  is  more 
limited.  Our  approach  is  specific  for  dimers;  examination  of  larger 
clusters  will  require  new  techniques.  Even  for  such  clusters,  however, 
a crude  estimate  of  the  distance  distribution  should  be  possible  by 
resorting  to  the  averaging  procedure  found  useful  for  dimers. 
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APPENDIX:  KOLMOGOROV'S  EQUATION  FOR  p 

X , 1 

If  the  primary  concern  is  with  the  probability  p^  that  site  x be 
occupied  regardless  of  cluster  configuration,  then  the  formalism  in 
Section  I suffices.  A more  elaborate  symbolism  is  required  to  evaluate 
the  probability  p . that  a specified  configuration  be  at  a given  site. 
In  general,  rate  constants  for  jumping  from  a site  to  its  neighbor  will 
depend  both  upon  the  initial  and  the  final  configuration  of  the  cluster. 

We  therefore  denote  the  rate  constant  for  transitions  from  x to  x+1 , 
starting  in  configuration  h and  ending  in  configuration  i.,  by 
similarly  the  rate  from  x to  x-1,  starting  in  configuration  h and  ending 
in  ^ is  Indicated  by  The  Kolmogorov  equation  for  the  probability 

p . of  having  a cluster  at  point  x in  configuration  i,  can  now  be  written 

X , 1 

as 


^Px.i 

dt 


“EX  i.,p  ..-E(X 

^ x-1  .hl'^x-l ,h  ' x,ih 

^ ^ P'x+l,hi  Px+l,h 


+ ^,ih^Px,i 


(Al) 


This  relation  is  crucial  to  the  discussion  of  trimer  motion,  for  example. 

It  is  the  basis  of  Eq.  (45),  for  the  probability  of  finding  a trimer  in 
configuration  lA. 

In  general,  this  detail  is  not  necessary.  We  then  regain  a relation 
for  the  probability  of  site  x being  occupied  , regardless  of  the  particular 
cluster  configuration,  by  summing  Eq.  (Al)  over  all  configurations  i at  x. 


That  is, 


dEp 


x,i 


x-1, hi  x-1, i ^ x,lh  x,lh  ‘^x,! 


(A2) 


dt 


^ P’x+l,hlPx+l,h 

h,l 
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The  general  approach  followed  here  is  similar  to  that  used  in  the  original 
derivations  of  the  probability  density  for  a randomized  random  walk.  See 
reference  5,  and  also  W,  Feller,  SIAM  (Soc.  Ind,  Appl.  Math.)  J.  Appl. 
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In  the  terminology  of  stochastic  analysis,  this  problem  corresponds  to 
a pure  birth  process. 
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Fig.  I 


Fig.  2 


Fig.  3 


Fig.  4 


Fig.  5 


FIGURE  CAPTIONS 

Schematic  of  dimer  and  trlmer  configurations  in  one-dimensional 
diffusion  on  a lattice  represented  by  shaded  circles.  Ooen 
circles  indicate  clusters.  Surface  sites  in  adjacent  channels 
are  throughout  assumed  to  be  exactly  abreast  of  each  other. 

Spacing  in  close-packed  direction  is  Z. 

Rate  constants  for  interconversion  of  dimer  configurations  on  a 
one-dimensional  lattice.  Configurations  are  labelled  by  the 
separation  between  atoms,  in  units  of  I,  along  the  direction  of 
diffusion. 

Configurations  allowed  on  the  points  x of  the  center  of  mass  lattice 
of  a dimer.  Arrows  indicate  the  jumps  and  jump  rates  between 
configurations  on  neighboring  sites. 

Schematic  of  potential  energy  relations  for  a dissociating  dimer. 

Shading  indicates  separation  at  which  interactions  between 

adatoms  in  cluster  vanish.  AE  denotes  diffusion  barrier  for  a 

a 

single  atom. 

Mean-square  displacement  of  dissociating  dimers,  modeled  on  the 

3 

behavior  of  Re  on  W(211).  Parameters  entering  estimates  are 
ttg  - 2,2  X 10^^  exp(-17,500/RT),  - 2.5  X 10^^  exp( -18 , 300/RT) , 

at  - 2,95  X 10^^  exp(-19,800/RT),  P^  - 10.5  expC-E^/RT), 

Pj/Pj  - exp(-E2/RT),  - 2.95  X 10^^  exp[-(19,800  + cEjI/RT]. 

All  estimates  at  T ■ 330°K.  Different  degrees  of  dissociation 
are  obtained  by  varying  Ej  between  0 and  5000  cal /mole. 
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Fig.  6. 


Fig.  7. 


Fig.  8. 


Fig.  9. 


Fig.  10. 


Fig.  11. 


Configurations  and  rate  constants  for  trimer  diffusion. 
Energetically  different  configurations  on  the  same  site  x are 
distinguished  by  A and  B. 

Probability  density  for  the  total  number  of  jumps  N made  by  a 
dimer  during  the  interval  t.  All  estimates  in  this  and  subsequent 
figures  are  made  for  the  same  effective  jump  rate  X - 2ab/(a+b), 
for  dimers  originating  in  even  positions. 

Probability  density  for  the  number  of  jumps,  averaged  over  adjacent 
values  of  N. 

t h 

Probability  density  for  the  time  at  which  the  jump  of  the 

dimer  occurs. 

Probability  density  for  the  displacement  x of  the  center  of 
mass  for  a dimer.  Continuous  curve  shows  probability  density 
for  Gaussian  walk  wit!i  same  effective  jump  rate. 

Probability  density  for  x,  averaged  over  adjacent  values  of  x. 
Density  for  Gaussian  walk  is  indicated  by  continuous  curve. 
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Fig.  1.  Schematic  of  dimer  and  trlmer  configurations  In  one-dlmenslonal 
diffusion  on  a lattice  represented  by  shaded  circles.  Open 
circles  Indicate  clusters.  Surface  sites  In  adjacent  channels 
are  throughout  assumed  to  be  exactly  abreast  of  each  other. 
Spacing  In  close-packed  direction  Is  I, 
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Fig.  2.  Rate  Constanta  for  Interconversion  of  dimer  configurations  on  a 
one-dimensional  lattice.  Configurations  are  labelled  by  the 
separation  between  atoms,  in  units  of  I,  along  the  direction  of 
diffusion. 
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Fig.  3.  Configurations  allowed  on  the  points  x of  the  center  of  mass  lattice 
of  a dimer.  Arrows  indicate  the  Jumps  and  Jump  rates  between 
configurations  on  neighboring  sites. 


Diffusion 


Interatomic  Potential 


Interatomic  Separation  (i) 


Schematic  of  potential  energy  relations  for  a dissociating  dimer, 
Shading  Indicates  separation  at  which  interactions  between 
adatoms  In  cluster  vanish.  denotes  diffusion  barrier  for  a 
single  atom. 
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Fig.  5.  Mean-square  displacement  of  dissociating  dimers,  modeled  on  the 

3 

behavior  of  Re  on  W(211).  Parameters  entering  estimates  are 


Oq  - 2.2  X 10 


12 

12 


exp(-17,500/RT),  p - 2.5  x 10^^  exp(-18,300/RT), 


1 


exp(-19,800/RT),  P - 10.5  exp(-E„/RT), 
®12  ^ 


a - 2.95  X 10 

P2/P1  - exp(-E2/RT),  - 2.95  x 10*-‘  exp[-(19,800  + cEjI/RT]. 
All  estimates  at  T - 330°K.  Different  degrees  of  dissociation 
are  obtained  by  varying  E2  between  o and  5000  cal/mole. 
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Fig.  6.  Configurations  and  rate  constants  for  trlmer  diffusion. 

Energetically  different  configurations  on  the  same  site  x ce 
distinguished  by  A and  B. 


Number  of  Jumps,  N 

Probabtltty  density  for  the  total  number  of 
dimer  during  the  interval  t.  All  estimates 
figures  are  made  for  the  same  effective  Jum| 
dimers  originating  in  even  positions. 


Gaussian 


Distance,  x • P 

Probability  density  for  the  displacement  x of  the  center  of  mass 
for  a dimer.  Continuous  curve  shows  probability  density  for 
Gaussian  walk  with  same  effective  jump  rate. 


